Cloning and expression profile of the alanine aminotransferase gene from kuruma shrimp Penaeus japonicus exposed to different salinities.

Several crustaceans including shrimps change the amount of specific free amino acids to regulate the osmotic pressure in their bodies. Kuruma shrimp Penaeus japonicus also increases the concentration of alanine (Ala) in the abdominal muscle following the increase of environmental salinity. In the present study, to elucidate the mechanisms of changes in Ala accumulation of kuruma shrimp depending on salinity, we cloned the gene encoding alanine aminotransferase (ALT), an enzyme involved in Ala biosynthesis, and examined its expression profile. It was found that the full-length kuruma shrimp ALT1 cDNA consisted of 3,301 bp, encoding 514 amino acids, and that all amino acid residues important for ALT activity were conserved. Phylogenetic analysis also indicated that the ALT gene cloned in this study was classified as ALT1. Moreover, we examined the expression levels of the ALT1 gene in the abdominal muscle and the hepatopancreas of kuruma shrimp acclimated at 1.7‰, 3.4‰, and 4.0‰ salinities, resulting that the mRNA levels of the ALT1 genes in both tissues of the shrimp acclimated at 4.0‰ were significantly higher than those at 1.7‰ for 12 h (p < 0.05). The mRNA levels of the ALT1 gene in the abdominal muscle of the shrimp acclimated for more than 24 h tended to increase following the increase of environmental salinity. These results indicate that ALT1 is responsible for the increase of free Ala concentration in the abdominal muscle of kuruma shrimp to regulate osmotic pressure at high salinity.

Metagenomic analysis of the microbial communities and associated network of nitrogen metabolism genes in the Ryukyu limestone aquifer.

While microbial biogeochemical activities such as those involving denitrification and sulfate reduction have been considered to play important roles in material cycling in various aquatic ecosystems, our current understanding of the microbial community in groundwater ecosystems is remarkably insufficient. To assess the groundwater in the Ryukyu limestone aquifer of Okinawa Island, which is located in the southernmost region of Japan, we performed metagenomic analysis on the microbial communities at the three sites and screened for functional genes associated with nitrogen metabolism. 16S rRNA amplicon analysis showed that bacteria accounted for 94-98% of the microbial communities, which included archaea at all three sites. The bacterial communities associated with nitrogen metabolism shifted by month at each site, indicating that this metabolism was accomplished by the bacterial community as a whole. Interestingly, site 3 contained much higher levels of the denitrification genes such as narG and napA than the other two sites. This site was thought to have undergone denitrification that was driven by high quantities of dissolved organic carbon (DOC). In contrast, site 2 was characterized by a high nitrate-nitrogen (NO3-N) content and a low amount of DOC, and this site yielded a moderate amount of denitrification genes. Site 1 showed markedly low amounts of all nitrogen metabolism genes. Overall, nitrogen metabolism in the Ryukyu limestone aquifer was found to change based on environmental factors.

Cold-induced muscle atrophy in zebrafish: Insights from swimming activity and gene expression analysis.

The investigation into the effects of cold acclimation on fish skeletal muscle function and its potential implications for muscle atrophy is of great interest to us. This study examines how rearing zebrafish at low temperatures affects their locomotor activity and the expression of genes associated with muscle atrophy. Zebrafish were exposed to temperatures ranging from 10 °C to 25 °C, and their swimming distance was measured. The expression levels of important muscle atrophy genes, Atrogin-1 and MuRF1, were also evaluated. Our findings show that swimming activity significantly decreases when the water temperature ranges from 10 °C to 15 °C, indicating a decrease in voluntary movement. Additionally, gene expression analysis shows a significant increase in the expression of Atrogin-1 and MuRF1 at 10 °C. This up-regulation could lead to muscle atrophy caused by decreased activity in cold temperatures. To investigate the effects of exercise on reducing muscle atrophy, we subjected zebrafish to forced swimming at a temperature of 8 °C for ten days. This treatment significantly reduced the expression of Atrogin-1 and MuRF1, emphasizing the importance of muscle stimulation in preventing muscle atrophy in zebrafish. These findings suggest that zebrafish can serve as a valuable model organism for studying muscle atrophy and can be utilized in drug screening for muscle atrophy-related disorders. Cold-reared zebrafish provide a practical and ethical approach to inducing disuse muscle atrophy, providing valuable insights into potential therapeutic strategies for addressing skeletal muscle atrophy.

Characterization of recombinant photoconverting green fluorescent Akanes.

Akanes are fluorescent proteins that have several fluorescence maxima. In this report, Akane1 and Akane3 from Scleronephthya gracillima were selected, successfully overexpressed in Escherichia coli and purified by affinity chromatography. Fluorescence spectra of the recombinant Akanes matured in darkness, or ambient light were found to have several fluorescence peaks. SDS-PAGE analysis revealed that Akanes matured in ambient light have two fragments. MS/MS analysis of Akanes digested with trypsin showed that the cleavage site is the same as observed for the photoconvertible fluorescent protein Kaede. The differences between the calculated masses from the amino acid sequence of Akane1 and the measured masses of Akane1 fragments obtained under ambient light coincided with those of Kaede. In contrast, a mass difference between the measured N-terminal Akane3 fragment and the calculated mass indicated that Akane3 is modified in the N-terminal region. These results indicate that numerous peaks in the fluorescent spectra of Akanes partly arise from isoproteins of Akanes and photoconversion. Photoconversion of Akane1 caused a fluorescence change from green to red, which was also observed for Akane3; however, the fluorescent intensity decreased dramatically when compared with that of Akane3.

Pilot study of a comprehensive resource estimation method from environmental DNA using universal D-loop amplification primers.

Many studies have investigated the ability of environmental DNA (eDNA) to identify the species. However, when individual species are to be identified, accurate estimation of their abundance using traditional eDNA analyses is still difficult. We previously developed a novel analytical method called HaCeD-Seq (haplotype count from eDNA by sequencing), which focuses on the mitochondrial D-loop sequence for eels and tuna. In this study, universal D-loop primers were designed to enable the comprehensive detection of multiple fish species by a single sequence. To sequence the full-length D-loop with high accuracy, we performed nanopore sequencing with unique molecular identifiers (UMI). In addition, to determine the D-loop reference sequence, whole genome sequencing was performed with thin coverage, and complete mitochondrial genomes were determined. We developed a UMI-based Nanopore D-loop sequencing analysis pipeline and released it as open-source software. We detected 5 out of 15 species (33%) and 10 haplotypes out of 35 individuals (29%) among the detected species. This study demonstrates the possibility of comprehensively obtaining information related to population size from eDNA. In the future, this method can be used to improve the accuracy of fish resource estimation, which is currently highly dependent on fishing catches.

Isolation and Characterization of Collagen and Collagen Peptides with Hyaluronidase Inhibition Activity Derived from the Skin of Marlin (Istiophoridae).

Type I and V collagens are the major components of fibrillogenic proteins in fish skin, and their hydrolysis products possess hyaluronidase inhibitory activity. In this study, for the first time, type I and V collagens were isolated from the skin of shortbill spearfish and striped marlin. Type I (2α1[I]α2[I]) and type V (α1[V]α3[V]α2[V]) collagens composed of distinct α-peptide chains with comparable structures were investigated using sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and UV spectrophotometric chromatography. After enzymatic digestion, the collagen peptides were purified by using ultrafiltration (30 KDa) and high-performance liquid chromatography (RP-HPLC) to yield CPI-F3 and CPV-F4 fractions with strong hyaluronidase inhibition rates (42.17% and 30.09%, respectively). Based on the results of simulated gastrointestinal fluid, temperature, and pH stability assays, CPI-F3 and CPV-F4 exhibited stability in gastric fluid and showed no significant changes under the temperature range from 50 to 70 °C (p > 0.05). The results of this first research on the bioactivity of type V collagen peptides provide valuable information for the biomedical industry and show the potential for future bioactivity investigations of type V collagen and its peptides.

A high-quality, haplotype-phased genome reconstruction reveals unexpected haplotype diversity in a pearl oyster.

Homologous chromosomes in the diploid genome are thought to contain equivalent genetic information, but this common concept has not been fully verified in animal genomes with high heterozygosity. Here we report a near-complete, haplotype-phased, genome assembly of the pearl oyster, Pinctada fucata, using hi-fidelity (HiFi) long reads and chromosome conformation capture data. This assembly includes 14 pairs of long scaffolds (>38 Mb) corresponding to chromosomes (2n = 28). The accuracy of the assembly, as measured by an analysis of k-mers, is estimated to be 99.99997%. Moreover, the haplotypes contain 95.2% and 95.9%, respectively, complete and single-copy BUSCO genes, demonstrating the high quality of the assembly. Transposons comprise 53.3% of the assembly and are a major contributor to structural variations. Despite overall collinearity between haplotypes, one of the chromosomal scaffolds contains megabase-scale non-syntenic regions, which necessarily have never been detected and resolved in conventional haplotype-merged assemblies. These regions encode expanded gene families of NACHT, DZIP3/hRUL138-like HEPN, and immunoglobulin domains, multiplying the immunity gene repertoire, which we hypothesize is important for the innate immune capability of pearl oysters. The pearl oyster genome provides insight into remarkable haplotype diversity in animals.

Environmental DNA study on aquatic ecosystem monitoring and management: Recent advances and prospects.

Environmental DNA (eDNA) is organismal DNA that can be detected in the environment and is derived from cellular material of organisms shed into aquatic or terrestrial environments. It can be sampled and monitored using molecular methods, which is important for the early detection of invasive and native species as well as the discovery of rare and cryptic species. While few reviews have summarized the latest findings on eDNA for most aquatic animal categories in the aquatic ecosystem, especially for aquatic eDNA processing and application. In the present review, we first performed a bibliometric network analysis of eDNA studies on aquatic animals. Subsequently, we summarized the abiotic and biotic factors affecting aquatic eDNA occurrence. We also systematically discussed the relevant experiments and analyses of aquatic eDNA from various aquatic organisms, including fish, molluscans, crustaceans, amphibians, and reptiles. Subsequently, we discussed the major achievements of eDNA application in studies on the aquatic ecosystem and environment. The application of eDNA will provide an entirely new paradigm for biodiversity conservation, environment monitoring, and aquatic species management at a global scale.

Age-Associated Different Transcriptome Profiling in Zebrafish and Rats: an Insight into the Diversity of Vertebrate Aging.

Most mammals, including humans, show obvious aging phenotypes, for example, loss of tissue plasticity and sarcopenia. In this regard, fish can be attractive models to study senescence because of their unique aging characteristics. The lifespan of fish varies widely, and several species can live for over 200 years. Moreover, some fish show anti-aging features and indeterminate growth throughout their life. Therefore, exploring the aging mechanism in fish could provide new insights into vertebrate aging. To this end, we conducted RNA sequencing (RNA-seq) assays for various organs and growth stages of zebrafish and compared the data with previously published RNA-seq data of rats. Age-associated differentially expressed genes (DEGs) for all zebrafish tissue samples reveal the upregulation of circadian genes and downregulation of hmgb3a. On one hand, a comparative analysis of DEG profiles associated with aging between zebrafish and rats identifies upregulation of circadian genes and downregulation of collagen genes as conserved transcriptome changes. On the other hand, in zebrafish, upregulation of autophagy-related genes in muscles and AP-1 transcription factor genes in various tissues is observed, which may imply fish-specific anti-aging characteristics. Consistent with our knowledge of mammalian aging, DEG profiles related to tissue senescence are observed in rats. We also detect age-associated downregulation of muscle homeostasis and differentiation-related genes in zebrafish gills, indicating a fish-specific senescence phenotype. Our results indicate both common and different aging profiles between fish and mammals, which could be used for future translational research.

Zebrafish Danio rerio myotomal muscle structure and growth from a spatial transcriptomics perspective.

Fish exhibit different muscle structures and growth characteristics compared with mammals. We used a spatial transcriptomics approach and examined myotomal muscle sections from zebrafish. Adult muscles were divided into eight regions according to spatial gene expression characteristics. Slow muscle was located in the wedge-shaped region near the lateral line and at the base of the dorsal fin, intermediate muscle was located in a ribbon-shaped region adjacent to slow muscle, and fast muscle was located in the deep region of the trunk, surrounded by intermediate muscle; the interior of fast muscle was further divided into 6 parts by their transcriptomic features. Combined analysis of adult and larval data revealed that adult muscles contain specific regions similar to larval muscles. These regions showed active myogenesis and a high expression of genes associated with muscle hyperplasia. This is the first study to apply spatial transcriptomics to fish myotomal muscle structure and growth.

Induction of endoplasmic reticulum stress markers in an acromegaly model.

Acromegaly is a growth hormone (GH) excess pathological condition in humans. Acromegaly is associated with somatic disfigurement and a wide range of systemic manifestations such as arthritis, neuropathy, carpal tunnel syndrome, reproductive disorders, metabolic disorders, and gastrointestinal complications. The influence of excess GH on the cellular level could aid in understanding the root causes of acromegaly-related health complications. Previously, we found that GH excess induces DNA damage to somatic cells and reduces the stem cells number and causes premature aging. In this study, an in-depth analysis of the acromegaly RNAseq data revealed the disruption of important biological cellular processes. Gene set enrichment analysis, heatmap, and enrichment analysis of acromegaly RNAseq data revealed induction of endoplasmic reticulum (ER) stress markers in various organs. Interestingly, the induction of ER stress was even more apparent than in aged zebrafish. Splicing of box-binding protein-1 (XBP1) mRNA is a hallmark of ER stress. Therefore, we quantified spliced XBP1 mRNA in different organs of our acromegaly model. Thus, our study emphasizes the importance of ER stress in GH oversecretion, which is important for understanding the health complications of acromegaly.

Exosome-derived small non-coding RNAs reveal immune response upon grafting transplantation in Pinctada fucata (Mollusca).

Exosomes, a subset of small extracellular vesicles, carry various nucleic acids, proteins, lipids, amino acids and metabolites. They function as a mode of intercellular communication and molecular transfer. Exosome cargo molecules, including small non-coding RNAs (sncRNAs), are involved in the immune response in various organisms. However, the role of exosome-derived sncRNAs in immune responses in molluscs remains unclear. Here, we aimed to reveal the sncRNAs involved in the immune response during grafting transplantation by the pearl oyster Pinctada fucata. Exosomes were successfully extracted from the P. fucata haemolymph during graft transplantation. Abundant microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs) were simultaneously discovered in P. fucata exosomes by small RNA sequencing. The expression patterns of the miRNAs and piRNAs at the grafting and initial stages were not substantially different, but varied significantly between the initial and later stages. Target prediction and functional analysis indicate that these miRNAs and piRNAs are related to immune response upon grafting transplantation, whereas piRNAs may also be associated with transposon silencing by targeting with genome transposon elements. This work provides the basis for a functional understanding of exosome-derived sncRNAs and helps to gain further insight into the PIWI/piRNA pathway function outside of germline cells in molluscs.

Visualisation of Phosphate in Subcalicoblastic Extracellular Calcifying Medium and on a Skeleton of Coral by Using a Novel Probe, Fluorescein-4-Isothiocyanate-Labelled Alendronic Acid.

The overload of nutrients of anthropogenic origin, including phosphate, onto coastal waters has been reported to have detrimental effects on corals. However, to the best of our knowledge, the phosphate concentration threshold for inhibiting coral calcification is unclear owing to a lack of information on the molecular mechanisms involved in the inhibitory effect of phosphate. Therefore, in this study, we prepared a new phosphate analogue, fluorescein-4-isothiocyanate (FITC)-labelled alendronic acid (FITC-AA), from commercially available reagents and used it as a novel probe to demonstrate its transfer pathway from ambient seawater into Acropora digitifera. When the juveniles at 1 d post-settlement were treated with FITC-AA in a laboratory tank, this phosphate analogue was found in the subcalicoblastic extracellular calcifying medium (SCM) and was absorbed on the basal plate in the juveniles within a few minutes. When the juveniles bear zooxanthellae at 3 months post-settlement, FITC-AA was observed on the corallite walls within a few minutes after adding ambient seawater. We concluded that FITC-AA in ambient seawater was transferred via a paracellular pathway to SCM and then absorbed on the coral CaCO3 skeletons because FITC-AA with a high polarity group cannot permeate through cell membranes.

Metagenomic analysis provides functional insights into seasonal change of a non-cyanobacterial prokaryotic community in temperate coastal waters.

The taxonomic compositions of marine prokaryotic communities are known to follow seasonal cycles, but functional metagenomic insights into this seasonality is still limited. We analyzed a total of 22 metagenomes collected at 11 time points over a 14-month period from two sites in Sendai Bay, Japan to obtain seasonal snapshots of predicted functional profiles of the non-cyanobacterial prokaryotic community. Along with taxonomic composition, functional gene composition varied seasonally and was related to chlorophyll a concentration, water temperature, and salinity. Spring phytoplankton bloom stimulated increased abundances of putative genes that encode enzymes in amino acid metabolism pathways. Several groups of functional genes, including those related to signal transduction and cellular communication, increased in abundance during the mid- to post-bloom period, which seemed to be associated with a particle-attached lifestyle. Alternatively, genes in carbon metabolism pathways were generally more abundant in the low chlorophyll a period than the bloom period. These results indicate that changes in trophic condition associated with seasonal phytoplankton succession altered the community function of prokaryotes. Our findings on seasonal changes of predicted function provide fundamental information for future research on the mechanisms that shape marine microbial communities.

Regulation of the Expression of the Myosin Heavy Chain (MYH) Gene myh14 in Zebrafish Development.

The human sarcomeric myosin heavy chain gene MYH14 contains an intronic microRNA, miR-499. Our previous studies demonstrated divergent genomic organization and expression patterns of myh14/miR-499 among teleosts; however, the regulatory mechanism is partly known. In this study, we report the regulation of myh14 expression in zebrafish, Danio rerio. Zebrafish myh14 has three paralogs, myh14-1, myh14-2, and myh14-3. Detailed promoter analysis suggested that a 5710-bp 5'-flanking region of myh14-1 and a 5641-bp region of myh14-3 contain a necessary regulatory region to recapitulate specific expression during embryonic development. The 5'-flanking region of zebrafish myh14-1 and its torafugu ortholog shared two distal and a single proximal conserved region. The two distal conserved regions had no effect on zebrafish myh14-1 expression, in contrast to torafugu expression, suggesting an alternative regulatory mechanism among the myh14 orthologs. Comparison among the 5'-flanking regions of the myh14 paralogs revealed two conserved regions. Deletion of these conserved regions significantly reduced the promoter activity of myh14-3 but had no effect on myh14-1, indicating different cis-regulatory mechanisms of myh14 paralogs. Loss of function of miR-499 resulted in a marked reduction in slow muscle fibers in embryonic development. Our study identified different cis-regulatory mechanisms controlling the expression of myh14/miR-499 and an indispensable role of miR-499 in muscle fiber-type specification in zebrafish.

Seasonal and annual changes in the microbial communities of Ofunato Bay, Japan, based on metagenomics.

Five years of datasets from 2015 to 2019 of whole genome shotgun sequencing for cells trapped on 0.2-µm filters of seawater collected monthly from Ofunato Bay, an enclosed bay in Japan, were analysed, which included the 2015 data that we had reported previously. Nucleotide sequences were determined for extracted DNA from three locations for both the upper (1 m) and deeper (8 or 10 m) depths. The biotic communities analysed at the domain level comprised bacteria, eukaryotes, archaea and viruses. The relative abundance of bacteria was over 60% in most months for the five years. The relative abundance of the SAR86 cluster was highest in the bacterial group, followed by Candidatus Pelagibacter and Planktomarina. The relative abundance of Ca. Pelagibacter showed no relationship with environmental factors, and those of SAR86 and Planktomarina showed positive correlations with salinity and dissolved oxygen, respectively. The bacterial community diversity showed seasonal changes, with high diversity around September and low diversity around January for all five years. Nonmetric multidimensional scaling analysis also revealed that the bacterial communities in the bay were grouped in a season-dependent manner and linked with environmental variables such as seawater temperature, salinity and dissolved oxygen.

Development of a time-series shotgun metagenomics database for monitoring microbial communities at the Pacific coast of Japan.

Although numerous metagenome, amplicon sequencing-based studies have been conducted to date to characterize marine microbial communities, relatively few have employed full metagenome shotgun sequencing to obtain a broader picture of the functional features of these marine microbial communities. Moreover, most of these studies only performed sporadic sampling, which is insufficient to understand an ecosystem comprehensively. In this study, we regularly conducted seawater sampling along the northeastern Pacific coast of Japan between March 2012 and May 2016. We collected 213 seawater samples and prepared size-based fractions to generate 454 subsets of samples for shotgun metagenome sequencing and analysis. We also determined the sequences of 16S rRNA (n = 111) and 18S rRNA (n = 47) gene amplicons from smaller sample subsets. We thereafter developed the Ocean Monitoring Database for time-series metagenomic data ( http://marine-meta.healthscience.sci.waseda.ac.jp/omd/ ), which provides a three-dimensional bird's-eye view of the data. This database includes results of digital DNA chip analysis, a novel method for estimating ocean characteristics such as water temperature from metagenomic data. Furthermore, we developed a novel classification method that includes more information about viruses than that acquired using BLAST. We further report the discovery of a large number of previously overlooked (TAG)n repeat sequences in the genomes of marine microbes. We predict that the availability of this time-series database will lead to major discoveries in marine microbiome research.

Phosphate bound to calcareous sediments hampers skeletal development of juvenile coral.

To test the hypothesis that terrestrial runoff affects the functions of calcareous sediments in coral reefs and hampers the development of corals, we analysed calcareous sediments with different levels of bound phosphate, collected from reef areas of Okinawajima, Japan. We confirmed that phosphate bound to calcareous sediments was readily released into ambient seawater, resulting in much higher concentrations of phosphorous in seawater from heavily polluted areas (4.3-19.0 µM as compared with less than 0.096 µM in natural ambient seawater). Additionally, we examined the effect of phosphate released from calcareous sediments on the development of Acropora digitifera coral juveniles. We found that high phosphate concentrations in seawater clearly inhibit the skeletal formation of coral juveniles. Our results demonstrate that calcareous sediments in reef areas play a crucial role in mediating the impact of terrestrial runoff on corals by storing and releasing phosphate in seawater.

N-Acetyl-d-Glucosamine-Binding Lectin in Acropora tenuis Attracts Specific Symbiodiniaceae Cell Culture Strains.

Many corals establish symbiosis with Symbiodiniaceae cells from surrounding environments, but very few Symbiodiniaceae cells exist in the water column. Given that the N-acetyl-d-glucosamine-binding lectin ActL attracts Symbiodiniaceae cells, we hypothesized that corals must attract Symbiodiniaceae cells using ActL to acquire them. Anti-ActL antibody inhibited acquisition of Symbiodiniaceae cells, and rearing seawater for juvenile Acropora tenuis contained ActL, suggesting that juvenile A. tenuis discharge ActL to attract these cells. Among eight Symbiodiniaceae cultured strains, ActL attracted NBRC102920 (Symbiodinium tridacnidorum) most strongly followed by CS-161 (Symbiodinium tridacnidorum), CCMP2556 (Durusdinium trenchii), and CCMP1633 (Breviolum sp.); however, it did not attract GTP-A6-Sy (Symbiodinium natans), CCMP421 (Effrenium voratum), FKM0207 (Fugacium sp.), and CS-156 (Fugacium sp.). Juvenile polyps of A. tenuis acquired limited Symbiodiniaceae cell strains, and the number of acquired Symbiodiniaceae cells in a polyp also differed from each other. The number of Symbiodiniaceae cells acquired by juvenile polyps of A. tenuis was correlated with the ActL chemotactic activity. Thus, ActL could be used to attract select Symbiodiniaceae cells and help Symbiodiniaceae cell acquisition in juvenile polyps of A. tenuis, facilitating establishment of symbiosis between A. tenuis and Symbiodiniaceae cells.

Estimation of tuna population by the improved analytical pipeline of unique molecular identifier-assisted HaCeD-Seq (haplotype count from eDNA).

Many studies have investigated the ability to identify species from environmental DNA (eDNA). However, even when individual species are identified, the accurate estimation of their abundances by traditional eDNA analyses has been still difficult. We previously developed a novel analytical method called HaCeD-Seq (Haplotype Count from eDNA), which focuses on the mitochondrial D-loop sequence. The D-loop is a rapidly evolving sequence and has been used to estimate the abundance of eel species in breeding water. In the current study, we have further improved this method by applying unique molecular identifier (UMI) tags, which eliminate the PCR and sequencing errors and extend the detection range by an order of magnitude. Based on this improved HaCeD-Seq pipeline, we computed the abundance of Pacific bluefin tuna (Thunnus orientalis) in aquarium tanks at the Tokyo Sea Life Park (Kasai, Tokyo, Japan). This tuna species is commercially important but is at high risk of resource depletion. With the developed UMI tag method, 90 out of 96 haplotypes (94%) were successfully detected from Pacific bluefin tuna eDNA. By contrast, only 29 out of 96 haplotypes (30%) were detected when UMI tags were not used. Our findings indicate the potential for conducting non-invasive fish stock surveys by sampling eDNA.